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1 Abstract 

° . 

We review the correspondence between strongly coupled lattice multiflavor Schwinger models 
and SU (Af) antiferromagnetic chains. We show that finding the low lying states of the gauge models 
is equivalent to solving an SU(AT) Heisenberg antiferromagnetic chain. For the two-flavor lattice 
Schwinger model the massless excitations correspond to gapless states of the Heisenberg chain, 
while the massive states are created by fermion transport in the ground state of the spin chain. 
Our analysis shows explicitly how spinons may arise in lattice gauge theories. 
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1 Introduction 



Though it is far from the scaling regime, the strong coupling limit is often used to study the 
qualitative properties of gauge field theories . Two important features of the spectrum of non- 
Abelian gauge theories appear there. The strong coupling limit exhibits confinement in a rather 
natural way. Furthermore with certain numbers of flavors and colors of dynamical quarks, it is 
straightforward to show that they also exhibit chiral symmetry breaking. This has motivated several 
more quantitative investigations of gauge theories using strong coupling techniques and there have 
been attempts to compute the mass spectrum of realistic models such as quantum chromodynamics 
(QCD) @ §. 

It has been recognized for some time that the strong coupling limit of lattice gauge theory 
with dynamical fermions is related to certain quantum spin systems. This is particularly true in 
the hamiltonian formalism and had already appeared in some of the earliest analyses of chiral 
symmetry breaking in the strong coupling limit Q . On the other side, it has been noted that there 
are several similarities between some condensed matter systems with lattice fermions and lattice 
gauge theory systems, particularly in their strong coupling limit. For example, it is well known 
that the quantum spin-1/2 Heisenberg antifcrromagnet is equivalent to the strong coupling limit of 
either U(l) or SU{2) lattice gauge theory ||, |, @, §, §. 

1 Lectures delivered by P. Sodano at the International Conference on "Mathematical Methods in Modern Theo- 
retical Physics" Tbilisi, Georgia, September 1998. 



Moreover, the staggered fermions which are used to put the Dirac operator on a lattice resemble 
ordinary lattice fermions used in tight binding models in condensed matter physics when the latter 
have an half-filled band and are placed in a background U(l) magnetic field it (mod 2ir) (1/2 of a 
flux quantum) through every plaquette of the lattice. This is an old result for d = 2 + 1 [[To] where 
it was already recognized in the first work on the Azbel-Wannier-Hofstaeder problem and has since 
been discussed in the context of the so-called flux phases of the Hubbard model. It is actually true 
for all d > 2 + 1 [Q, i.e. staggered fermions with the minimal number of flavors are identical to 
a simple nearest neighbor hopping problem with a background magnetic field which has 1/2 flux 
quanta per plaquette. In the condensed matter context, the magnetic flux can be produced by a 
condensate, as in the flux phases. As an external magnetic field, for ordinary lattice spacing it is as 
yet an experimentally inaccessible flux density. However, it could be achieved in analog experiments 
where macroscopic arrays of Josephson junctions, for example, take the place of atoms at lattice 
sites. 

In (l2) it was shown that the strong coupling limit of d-dimensional QED with 2 d /2^ d / 2 ^ flavors of 
fermions can be mapped onto the S — 1/2 quantum Heiscnberg antiferromagnet in d— 1 dimensions. 
The chiral symmetry breaking phase transition corresponds to a transition between the flux phase 
and the conventional Neel ordered phase of an antiferromagnet. 

In this lectures we shall discuss the correspondence between lattice gauge theories and spin 
systems by analyzing, as a concrete example, the mapping of the strongly coupled multiflavor lattice 
Schwinger models — i.e. QED2 with massless fermions Jl3| — into the antiferromagnetic quantum 
spin-1/2 Heisenberg chains. In the one-flavor Schwinger model, the anomaly is realized in the lattice 
theory via spontaneous breaking of a residual chiral symmetry p[ . In the multiflavor models, at 
variance with the one-flavor case, the chiral symmetry is explicitly broken by the staggered fermions, 
and the non-zero vacuum expectation value of a pertinent condensate is the only relic on the lattice 
of the chiral anomaly in the continuum |l^, [ll| . 

In section 2 we introduce the antiferromagnetic quantum spin-1/2 chains; we review the Bcthc 
ansatz Jl7| solution of the Heisenberg chain jl8| and compare it with the exact diagonalization of 4, 6 
and 8 sites chains. Furthermore we generalize the analysis to the spin-1/2 SU(Af) antiferromagnetic 
Heisenberg chains |l9|, [l5|, [ll| . 

Section 3 is devoted to a detailed investigation of the two-flavor lattice Schwinger model in the 
strong coupling limit using the hamiltonian formalism and staggered fermions. We show that the 
problem of finding the low lying states is equivalent to solving the Heisenberg antiferromagnetic spin 
chain. The massless excitations of the gauge model correspond to gapless states of the Heisenberg 
chain, while the massive excitations are created by fermion transport in the ground state of the spin 
chain and we compute their masses in terms of vacuum expectation values (VEV's) of spin-spin 
correlation functions |§ ||, Q. 

In section 4 we study the generic A/"-flavor lattice Schwinger models and their mapping into 
spin-1/2 quantum SU(Af) antiferromagnetic spin chains. The analysis strictly parallels what we 
found in section 3, but some rather surprising difference arises depending on if Af is even or odd. 



2 Antiferromagnetic spin chains 

The title of this section deserves an entire book to be properly developed, and in fact some books, 
reviews and long articles about spin chains do already exist pp[ . Spin chains interested mathe- 
matical physicists for their exact solvability, field theorists for the possibility to test their methods 
by comparation with exact results and condensed matter physicists for the possibility to under- 
stand the spectra of real physical systems. It is our purpose to report a self-contained analysis of 
antiferromagnetic spin chains, useful to understand section 3 and 4 where we shall show that the 
multiflavor lattice Schwinger models in the strong coupling limit are effectively described by the 
spin-1/2 SU(AT) antiferromagnetic Heisenberg chains. 

A common tool to study spin systems is spin-wave theory. Spin-wave theory was developed 
by Anderson, Kubo and others [pl| in the 1950s and it is a very powerful approach to quantum 
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magnetic systems in dimensions greater than one, predicting long-range order and gapless Goldstonc 
bosons. The situation remained clouded for magnetic chains because it was known that long-range 
order could not occur in one dimension p2| , but the Bethe ansatz predicts massless excitations. 
There is "quasi-long-range-order" corresponding to a power law decay of spin-spin correlators and 
there are gapless excitations which are not true Goldstone bosons. 

Quantum spin chains have been extensively studied in the literature, starting from the seminal 
paper by Hans Bethe jl7j in 1931 for the spin-1/2 case, where he introduced an ansatz to write 
down the eigenfunctions of the Hciscnbcrg Hamiltonian, describing a chain with periodic boundary 
conditions. Seven years later Hulthen |2^] was able to compute the ground state energy of the 
antiferromagnetic Heisenberg chain. We had to wait until 1984 to know exactly the spectrum of 
this model, when Faddeev and Takhtadzhyan Jl8| analyzed the model using the algebraic Bethe 
ansatz and showed that the only one-particle excitation is a doublet of spin-1/2 quantum excitations 
with gapless dispersion relation. This excitation is a kink rather than an ordinary particle and is 
called spinon. All the eigenstates of the antiferromagnetic Heisenberg Hamiltonian contain an even 
number of kinks, nevertheless the kinks arc localizable objects and one can consider their scattering. 
There are no bound states of kinks in this model. 

The one dimensional spin systems are not only unusual because they are exactly solvable but 
also for their incompatibility with long range order. Actually solvability and absence of long range 
order are deeply related concepts; systems with spontaneously broken symmetries are more difficult 
to describe and resist analytical solutions. All isotropic antiferromagnetic quantum spin chains 
with short range interactions exhibit quantum spin-liquid ground states — i.e. states with short 
range antiferromagnetic correlations and no order. Moreover these chains have strange quantum 
elementary excitations above these ground states, that are not ordinary spin waves but are usually 
called spinons — neutral spin-1/2 kinks. 

As R. B. Laughlin J24|] points out: "spinons are not "like" anything familiar to most of us, but are 
instead an important and beautiful instance of fractional quantization, the physical phenomenon 
in which particles carrying pieces of a fundamental quantum number, such as charge or spin, arc 
created as a collective motion of many conventional particles obeying quantum mechanical laws. 
The fractional quantum number of the spinon is its spin. It is fractional because the particles out 
of which the magnetic states are constructed are spin flips, which carry spin 1." 

One significant result of our analysis Jl9|, [lj| has been to show explicitly how spinons appear 
in lattice gauge theories. More precisely we showed that the massless excitations of the two-flavor 
Schwinger model coincide with the spinons of the spin-1/2 antiferromagnetic Heisenberg model. 



Section (2.1) is devoted to review the Bethe ansatz solution of the antiferromagnetic Heisenberg 
chain. We shall compare the exact solution given in [[l8) with a study of finite size chains of 4,6, 
and 8 sites in section (2.2). An original result presented is the thermodynamic limit coefficient of 
the states with N — 2 domain walls composing the ground state. 

The spin-spin correlation function < g.s.\S x ■ S x +2\g-S- > originally computed by M. Taka- 
hashi |2f|, is discussed in section (2.3). The vacuum expectation value of the square of the vector 
V = J2 X A Sx+i is a l so computed. 

In section (2.4) we introduce the spin-1/2 SU(JV) antiferromagnetic Heisenberg chains. To 
provide an intuitive picture we study the SU(3) two sites chain and we find the ground state. A 
short review of the literature on SU(Af) antiferromagnetic Heisenberg chains is provided. 



2.1 The Bethe Ansatz solution of the antiferromagnetic Heisenberg chain 

In this section the Bethe ansatz solution of the antiferromagnetic Heisenberg chain jl8| is dis- 
cussed in detail. Moreover we shall compare the exact solution given by L.D. Faddeev and L.A. 
Takhtadzhyan in 111] with a study of finite size chains of 4, 6 and 8 sites. We show that already 
these very small finite systems exhibit spectra that match very well with the thermodynamic limit 
solution. We suggest to the reader interested in the subject the references [p0| . 
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The Bethe ansatz is a method of solution of a number of quantum field theory and statistical 
mechanics models in two space-time dimensions. This method was first suggested by Bethe |l7[] 
in 1931 from which takes its name. Historically one can call this formulation the Coordinate 
Bethe ansatz to distinguish it from the modern formulation known as Algebraic Bethe ansatz. 
The eigenfunctions of some (l+l)-dimensional Hamiltonians can be constructed imposing periodic 
boundary conditions which lead to a system of equations for the permitted values of momenta. 
These are known as the Bethe equations which are also useful in the thermodynamic limit. The 
energy of the ground state may be calculated in this limit and its excitations can be investigated. 



We review here the method applyied to the study of the antiferromagnetic Heisenberg chain 18 1. 
In particular it is shown that there is only one excitation with spin-1/2 which is a kink: physical 
states have only an even number of kinks, therefore they always have an integer spin. The one 
dimensional isotropic Heisenberg model describes a system of N interacting spin-^ particles. The 
Hamiltonian of the model is 

N 1 
Hj = jJ2(S x -S x+1 --) . (2.1) 

x=l 

where J > ( J < would describe a ferromagnet) and the spin operators have the following form 

S x = lr ® 1 2 <g> y <g> ljv . (2.2) 

They act nontrivially only on the Hilbert space of the X th site. Periodic boundary conditions are 
assumed. 



The Hamiltonian (ETj) is invariant under global rotations in the spin space, generated by 



S = Y,Sx ■ (2.3) 



x = l 

Due to the periodic boundary conditions, under translations generated by the operator T, 

f :S X ^ S x _i (2.4) 

the Hamiltonian is invariant and [S'jT] = 0. 

In order to diagonalize Hj it is convenient to use an eigenfunction basis of operators commuting 
with Hj, so obviously S 2 , S z and also T. Let us sketch the Coordinate Bethe ansatz technique. 
One has to introduce the "false vacuum" 

N 



ift>=n®iT>* ( 2 - 5 ) 



x=l 



with 



S+\ T>* = (2.6) 



5||T>» = JlT>x (2.7) 



where 



Si ± S 2 X (2.8) 



S 3 \fl> = y|^> (2.9) 

s 2 \n> = ^(| + i)|o> (2.io) 
f\n> = \n> . (2.ii) 
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All the other basis vectors have S 3 < ^ and one can get them by properly acting on |fl > with the 
lowering operators S~ . Let us start from the generic state with S 3 = f - 1, with N-l spins up 
and M = 1 spin down 



JV 



\M = 1 >= ^4> x \x > with \x>=S x \Q> 



(2.12) 



where the coefficients <ft x must be such that \M = 1 > is a translationally and rotationally invariant 
state 



A' 



N N 

f\M = 1 >=^2(j) x \x - 1 >=^<j) x+l \x >=nY^<j> x \x > 



and from Eq.( [2.13 ) one gets 



x=l 



x = l 



4>x+l = A*0x 



(2.13) 



(2.14) 



b x+x = fj, x (j)i x ^ N 
4>i = = ^ N( t>i — * ^ N = 1 

,x-l 



(2.15) 
(2.16) 



Setting <f>x = 1 one has (j> x = [i x . There are N possible values for 4> x , corresponding to the N 



N 



roots of the unity. One of these roots corresponds to a state with S — , 



N 



J2\ x >= S~\Q> — > n = l 



(2.17) 



x=l 



while the other N — 1 roots have S = y — 1. 

The generic case with M spins down is more complicated. Let us consider the case M = 2 to 
understand what happens 



|M = 2>= cj)(xi,X2)\xi,X2 > with \x\, X2 >— S Xl S X2 \£l > 

xi <x 2 = l 

By requiring the translational invariance of the state \M = 2 > one has 

f\M = 2 >= T\M = 2 > 

and for X2 < N one has 



(xi + 1,X2 + 1) = T(f)(xi,X2) 



(2.18) 

(2.19) 

(2.20) 
(2.21) 

but due to the periodic boundary conditions one has to find coefficients <j){x\,X2) that satisfy not 
only Eq. fl2.20D but also 

4>(l,x+ 1) = T<f>(x, N) . (2.22) 

Eq.(2.22) is no more satisfied by ( p. 21 ). Bethe proposed the following ansatz for the coefficients 

(j)(xi,x 2 ) 

cf>(x u x 2 )=A lr2 ^- 1 ^- 1 + A 2tl ^- 1 ^- 1 . (2.23) 
Eq.( 2.20| , 2.22 ) are satisfied if the following equations hold 



that would easily give 



zi,%2) = Mi A*2 > T = M1M2 



A 12 = A 2 itf 

A 2l = A 12 pL2 

By imposing to \M = 2 > to be an highest weight state 

S + \M = 2 >= 



(2.24) 
(2.25) 



(2.26) 
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taking into account Eq. ( |2.24 2.25Q and introducing the following change of variables 

\ _ l 



one gets the so called Bethe ansatz equations 



, A Q 



2 \N _ TT A a \p — i 



In the general case of M spins flipped one gets 



(rr*)"=-n 



Aq — A/3 — j 

A Q - A/3 + t 



(2.27) 



(2.28) 



(2.29) 



The energy and the momentum of a given state with M spins down can be expressed in terms 
of the parameters X a 



E 



M 



M 

E- 



2 ^ A^ 

a— 1 a 



i 



M 



M 



= i In T = ^2 Pa = i 



M 



(2.30) 



(2.31) 



Energy and momentum are thus additive as if there were M independent particles and the A Q 
must satisfy the Bethe ansatz equations (2.29) in order for Em and Pm to be eigenvalues of the 
Hamiltonian and momentum operators. 

The solution of the antiferromagnetic Heisenberg chain is reduced to the solution of the system 
of the M algebraic equations ( 2.29| ). This, in general, is not an easy task. It can be shown |18|, 
however, that, in the thermodynamic limit N — > oo, the complex parameters A have the form 



A; = Xj,L + H 



I 



-L,-L + l, 



1,L: 



(2.32) 



where L is a non- negative integer or half-integer, Xj.L is the real part of the solution of ( 2.29 ) and 
we shall define shortly the set of allowed values for the integer index j. The A's that, for a given 
\j,L, are obtained varying I between \—L,L] by integer steps, form a string of length 2L + 1, see 
fig. ([!]). This arrangement of A's in the complex plane is called the "string hypothesis" In 
the following we shall verify that, even on finite size systems, the "string hypothesis" is very well 
fulfilled. 

In a generic Bethe state with M spins down, there are M solutions to ( 2.29| ), which can be 
grouped according to the length of their strings. Let us denote by vl the number of strings of 
length 2L + 1, L = 0, i, . . .; strings of the same length are obtained by changing the real parts, A^, 
of the A's in (2.32); as a consequence j — 1, . . . , v^. If one denotes the total number of strings by 
q one has 

q = J2»L ' M = J2( 2L + 1 >l ■ (2.33) 



The set of integers (M, q, {vl}) constrained by (2.33), characterizes Bethe states up to the 
fixing of the q numbers \j,L', this set is called the "configuration". Varying M, q and Ul, one is 
able to construct all the 2 N eigenstates of an Heisenberg antiferromagnetic chain of N sites |18 . 
The energy and momentum of the Bcthe's state, corresponding to a given configuration — within 
exponential accuracy as N — * oo — consist of q summands representing the energy and momentum 
of separate strings. For the parameters A^l of the given configuration, taking the logarithm of 
(2.29) the following system of equations is obtained in the thermodynamic limit 



2N arctan ■ 



A 



In 



2-kQ 



vl. 



i,Li 



L 2 fe=l 



(A 



3,Li 



,) 



(2.34) 
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Figure 1: Strings for L = 0, 5, 1, | 



where 



L1+L2 , ^ 

< I ) l 1 l,(A) = 2 > (arctan — + arctan ) 

L L + l 

L=\L!~L 2 \^0 



(2.35) 



Integer and half integer numbers Qj.L parametrize the branches of the arcotangents and, con- 
sequently, the possible solutions of the system of Eqs.(2.34). In ref. [jl8| it was shown that the Qj.L 
arc limited as 



with Q™ ax given by 



QT X < Qi,l < Q2,l <■■■< Qv l ,l < QT 



Q 



rn a x 
L 



(2.36) 
(2.37) 



J(L U L 2 ) 



(2.38) 



and 

2mm(Li,L 2 ) + l if L x ^ L 2 
2Li + \ if L x = L 2 . 

The admissible values for the numbers Q^l are called the "vacancies" and their number for every 
L is denoted by Pl 

P L = 2Q n L iax + 1 . (2.39) 

The main hypothesis formulated in (l8) is that to every admissible collection of Qj.L there 
corresponds a unique solution of the system of equations ( [2.34 ). The solution always provides, in a 
multiplet, the state with the highest value of the third spin component S 3 . 

Let us now consider some simple example. The simplest configuration has only strings of length 
1, i.e. all the A's are real. The singlet associated to this configuration 



N 

M = q = is a = — , vl = 



L > 



(2.40) 



is the ground state. The vacancies of the strings of length 1, i.e. the admissible values of Qj.o, due 
to eqs.( p.36p37| , |2.38[ ), belong to the segment 



N 



1 



1 ^ N 
h - < Qa o < 

4 2 - 4 2 



(2.41) 
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Therefore they are N/2. All these vacancies must then be used to find the N/2 strings of length 1. 
As a consequence this state is uniquely specified and no degeneracy is possible. 

Next we consider the configuration that provides a singlet with 1 string of length 2 and all the 
others of length 1: 

N N N 1 

M = - , g = - - 1 , u = --2 , v, = 1 , v L = , L > - . (2.42) 

For the strings of length 1 the number of vacancies is again N/2; for the string of length 2 there is 
one vacancy and the only admissible Qj t i equals 0. Thus, since the number of strings of length 1 is 



^o = ^ — 2, there are two vacancies for which Eqs.(2.34) have no solution; they are called "holes" 



2 

and are denoted and Q^. This configuration is determined by two parameters: the positions 



of two "holes" which vary independently in the interval (2.41) 



There is another state with only 2 holes: the triplet corresponding to the configuration 

N 

M = q=v Q = — -l , v L = , L > (2.43) 
The number of vacancies for the strings of length 1 equals y + 1 , while v§ = ~ — 1 . 



The excitations determined by the configurations (2.42 2.43) belong to the configuration class 



called in |18j M.at- The class M.at is characterized as follows: the number of strings of length 
1 in each configuration belonging to this class differs by a finite quantity from N/2, Vq = y — fco 
where k is a positive finite constant, so that the number of strings of length greater than 1 is finite. 
From ( |2~39| ) one then has 



N k - 

P = —+ko-2^v L (2.44) 



2 

L>0 



P L = 2k -2 J{L,L')v L , , L>0 (2.45) 



L'>0 

so that 



N 

Po > y , Pl < 2k , L > . (2.46) 
From ( [2.44 ) follows that the number of holes for the strings of length 1 is always even and equals 



2 only for the singlet and the triplet excitations discussed above. One can imagine the class M.at 
as a sea of strings of length 1 with a finite number of strings of length greater than 1 immersed 
into it. It was proven in p"8| that, in the thermodynamic limit, the states belonging to Mat have 
finite energy and momentum with respect to the antiferromagnetic vacuum, whereas each of the 
states which corresponds to a configuration not included in the class M. at nas an infinite energy 
relative to the antiferromagnetic vacuum. 



Let us now sketch the computation of the thermodynamic limit ground state energy. Eqs.(2.34) 
for the ground state have the form 



■kQ- 1 N/2 

arctan2A :) = ^ + — arctan(Aj — A&) . (2.47) 



fc=i 

Taking the thermodynamic limit N — > oo, one has 



9l^x , -\<*<\ ■ A, -AI.D . (2.4N, 
and Eqs.( 2.47 ) can be rewritten in the form 

arctan2A(x) = irx + / arctan(A(cc) — X(y))dy . (2.49) 
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Upon introducing the density of the numbers A(x) in the interval dX 



= W (2 - 50) 

dx \x=x(\) 

and differentiating Eqs.( 2.49| ), one gets 

p(A ) = J- r _L^ e -^lde= -4— . (2.51) 
^ ; 27r./_ 00 1 + e-ISI 2cosh7rA v ; 

The density p(X) introduced in this way is normalized to 1/2. It is now easy to compute the energy 
and the momentum of the ground state 

N_ 

2 poo jrjy pOO j 

E g . s . = J2^= N l e(X)p(X)dX = -— I ^ (A2 + ^ )cosh7rA =-JiVln2 (2.52) 



Pg.s. = Y / Pa=N P (X)p(X)dX = - 1 r dX - = -7T (mod 2tt) . (2.53) 

~l J-oo 2 J -co coshTrA 2 

According to Eq.( 2.53| ), P g . s . = (mod27r) for ^ even, and P g . s . = it (mod 27r) for ^ odd. The 
ground state, as expected, is a singlet, in fact the spin S is given by 



N ^ N 



a=l 



S=--^l = --N / PWd^O • (2.54) 



Let us analyze the triplet described by Eq.( [2.43 ); Eqs.(2.34) take the form 

4-1 

arctan2A :) ' = — — - + — ^ arctan(Aj — Afc) (2.55) 
k=\ 

where now the numbers Qj lie in the segment [-O] and have two holes, Q± and with 
Qi < ■ Taking the thermodynamic limit one gets 

o {h) o {h) O 1 

, , ^^x+-(e(x-x 1 ) + ^- 2 ; 2 )) (2.56) 



where 6{x) is the Heaviside function. Eqs.( |2.55|) become 

arctan2A(x) = irx + -^{9(x — X\) + 8(x — £2)) + / arctan(A(x) — X(y))dy . (2.57) 



Eq. fl2.57D gives, for this triplet, the density of A, p(X) = |j 



p t (A) = p(X) + l(a(A - Ai) - a(A - A 2 )) (2.58) 
where p(A) is given in Eq.(pT5ll) and 

°w=- toLu^=® e « ■ (2 - 59) 

Ai and A 2 are the parameters of the holes, A^ = X(xi), i = 1,2. The energy and the momentum of 
this state measured from the ground state are now easily computed 

P 00 

e T (Ai, X 2 )=N e{X){p t {X) - p{X))dX = e(A x ) + e(A 2 ) (2.60) 
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p T (Xi,X 2 ) = N I p(X)(p t {X)-p(X))dX=p(X 1 )+p(X 2 ) (mod 2tt) (2.61) 

where 



<X) = J_J,MX-,W = J J -^ X (2.62) 

/oo 
7T 

p(p)a(X — n)dp — arctan sinh 7T A — — , — n < p(X) < . (2.63) 
-oo ^ 

From Eqs.( |2.62| , ^63|) one gets 

Jn 

e = — Y sin P ■ ( 2 - 64 ) 

The momentum pt(X±, X 2 ) varies over the interval [0, 2ir), when Ai and A2 run independently over 
the whole real axis. The spin of this state can be computed by the formula 

/oo 
(<t(A - Ai) + a(X - X 2 ))dX = 1 . (2.65) 
-00 

Let us finally analize the singlet excitation characterized by the configuration ( 2.42| ). Denoting 
by Ag the only number among the Xj^/2 which characterizes the string of length 2 and by Xj the 
numbers Xj.o for the strings of length 1, Eqs. fl2.34 ) read 



arctan 2Aj 



— —2 

7T<2j 



jU(Aj - A s ) + ^ £ arctan(A, - A fc ) (2.66) 



N N N 

fc=i 



with 



f -2 

arctanAs = — $(A S - Xj) (2.67) 

3=1 

2 

$(A) = arctan 2A + arctan -A . (2.68) 

3 

The ^ — 2 numbers Qj vary in the segment [— + 4, ^ — i]; among them there are the two holes 
Qi and Q 2 . Taking the thermodynamic limit one finds the density of A's for the singlet 

p(As) - p(X) + ^(ct(A - Ai) + <r(A - A 2 ) + c^(A - A s )) (2.69) 



where p and a were given in Eqs.(|||,[||l) and where 



1 f°° 9 

"(A) = -^/ e-*lt'-^dg = - ; ■ (2.70) 
2irJ_ QO tt(1 + 4A 2 ) 

In JlSj it was demonstrated that the string parameter As is uniquely determined by the A's 
parametrizing the two holes 

As= Al + Az ■ (2.71) 

In it was also proved the remarkable fact that the string of length 2 does not contribute to the 
energy and momentum of the excitation, so that the singlet and the triplet have the same dispersion 
relations 

es(Ai,A 2 ) = eT(Ai,A 2 ) = e(A 1 ) + e(A 2 ) (2.72) 

Ps(Xi,X 2 ) = p T (Xi,X 2 )=p(X 1 )+p(X 2 ) (mod 2tt) . (2.73) 

The spin of this excitation is, of course, zero 

poo 

S=-2- (2<r(A) +uj(X))d\ = (2.74) 
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The only difference between the state whose configuration is given in Eq.(2.43) and the state of 
Eq.( |2.42D is the spin. 

To summarize, the finite energy excitations of the antiferromagnetic Heisenberg chain are only 
those belonging to the class Ai af an d are described by scattering states of an even number of 
quasiparticlcs or kinks. The momentum p of these kinks runs over half the Brillouin zone — tt < p < 
0, the dispersion relation is e(p) = 4^ ship, Eq.(2.64), and the spin of a kink is 1/2. The singlet and 
the triplet excitations described above are the only states composed of two kinks, the spins of the 
kinks being parallel for the triplet and antiparallel for the singlet. For vanishing total momentum 
all the states belonging to Mat have the same energy of the ground state so that they are gapless 
excitations. Since the eigenstates of Hj always contain an even number of kinks, the dispersion 
relation is determined by a set of two-parameters: the momenta of the even number of kinks whose 
scattering provides the excitation. There are no bound states of kinks. 



2.2 Finite size antiferromagnetic Heisenberg chains 

Let us now turn to the computation of the spectrum of finite size quantum antiferromagnetic 
chains by exact diagonalization. We shall see that already for very small chains, the spectrum is 
well described by the Bethe ansatz solution in the thermodynamic limit. Furthermore, an intuitive 
picture of the ground state and of the lowest lying excitations of the strongly coupled two-flavor 
lattice Schwinger model emerges, due to the mapping of the gauge model onto the spin chain — see 
section 3. 

The states of an antiferromagnetic chain are classified according to the quantum numbers of 
spin, third spin component, energy and momentum \S, S 3 ,E,p >. For a 4 site chain the momenta 
allowed for the states are: 0, §, ^ mod 2ir. The ground state is 

\g. s . >= 10,0,-37,0 >= -L( 2 I UU> +2| iUT> -| TTU> -| TUT> -i un> -I ITU>) • 

V 12 



(2.75) 

This state is P-parity even. In fact, by the definition of P-parity given in Eq.( 3.178 ), the P-parity 
inverted state (2.7E) is obtained by reverting the order of the spins in each vector | . . . > appearing 

in (H), e.g. | UTT>^ I TTU>- 

The A's associated to the ground state (solution of the Bethe ansatz equations ( 2.29| )) arc 
At = =7= and A2 = —K=. There is also an excited singlet 

2v3 2v3 

|o,o,-J,7r>=4f(lUn>-Um>-|TUT>+ITTU>) ■ (2.76) 

V4 

It is P-even, so that it is a S p = + excitation, with the same quantum numbers (the isospin 
is replaced by the spin) of the lowest lying singlet excitation of the strongly coupled Schwinger 



model discussed by Coleman pM. The state (2.76) also coincides with the excited singlet described 



by the configuration (2.42). It has only two complex A's which arrange themselves in a string 



approximately of length 2, Ai = — A2 = iy 17 and there are two holes with = — | and 
n W _ 1 

W2 — 2" 

There are also three excited triplets, whose highest weight states are 

|i,i,-J,|> = ^(UtTt>+*|UTt>-|mt>-i|tm>) (2.77) 

|i,i,-2J,7r> = i=(|4tTT>-|TITT>+|TTIT>-|TTTI>) (2.78) 

V4 

iu,-j,y> = ^|(um>-*iTm>-imT>+*iTm>) • (2.79) 

Among these, only the non-degenerate state with the lowest energy has a well defined P-parity 



(2.78). It is a S p — 1 like the lowest lying triplet of the two-flavor strongly coupled Schwinger 



model. The degenerate states can be always combined to form a P-odd state. 
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We thus see that within the states in a given configuration there is always a representative state 
with well defined parity, the others are degenerate and can be used to construct state of well defined 
energy and parity. Moreover the parity of the representative states (with respect to the parity of 
the ground state) is the same of the one of the lowest-lying Schwinger model excitations in strong 
coupling. 

All the triplets in fl2.79| ) have one real A and two holes; they can be associated with the family 
of triplets ( f2~43| ). In table (Q) we summarize the triplet A's and Q^'s. The spectrum exhibits also 



Table 1: Triplet internal quantum numbers 



TRIPLET 


A 




|l,l,-J,f > 


1 

2 


-1 


|1,1,-2J,tt> 





-l i 


|M,-J,tF > 


1 

2 


1 



a quintet, whose highest weight state is 

12,2,0,0 >= | TTTT> 

We report in fig.(@) the spectrum of the 4 sites chain. 



E/J 



[3- 



tt/2 



3tt/2 2tt 



& B ^ 

© 



■a 



O QUINTET 
Q TRIPLETS 
□ SINGLETS 



P 



(2.80) 



Figure 2: Four sites chain spectrum 

Let us analize the spectrum of the 6 site antiferromagnetic chain. The momenta allowed for the 
states are now 0, |, ^r, n, ^ mod 27r. The ground state is 

!<?•*•> = |0,0,-^(5 + Vl3),7r>=-=i={|mUT>-|UTm> 
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+ 1 mui> -i utut> +i mm> -i mm> +i umt> -i mm> 
- i utm> +i mm> -i mtu> +i imiT> -i mm> +i umt>) 

+ 4 ~ 3 (i rmu> -i mur> +i uun> -i mm> +i unu> -i mm>)} (2.si) 

This state is odd under P-parity. The spectrum of the six sites chain is reported in fig.(^). There 
are 9 triplets in the spectrum. In it was already pointed out that the number of lowest lying 
triplets for a finite system with N sites is N(N + 2)/8, so for N = 6 there are 6 lowest lying triplet 
states. In order to identify these 6 states among the 9 that are exhibited by the spectrum of fig.(||), 
it is necessary to compute their A's and their Q's. In this way in fact, We can find o ut which are the 
triplets characterized by two holes and thus belonging to the triplet of type ( [2.43 ). In table (||) we 



report the internal quantum numbers of the lowest lying triplets. The Q^ h ''s vary in the segment 
[—§,§]• The highest weight state of the triplet of zero momentum and energy — (J/2)(5 + V5) reads 

io,o,-^(5 + %/5),o> - - 1 r 3 t^ dutm>+i mm>+i mtu>+i mut> 

2 V45 - 15\/5 2 

+ i mm> +i rum>) 

+ (i mm> +i ntm> +i mut> +i mm> +i mur> +i utm>) 

+ (l - V5)(| ITTITT> +1 TTITTI> +1 TITTIT>)} ■ ( 2 -82) 



One can get the triplet of energy — ( J/2)(5 — ^/5) from (2.82) by changing V5 — ► — \/5. As can be 



explicitly checked from (2.82), the two non-degenerate triplets of zero momentum are then P-parity 
even, namely they have opposite parity with respect to that of the ground state, as it happens 
for the lowest lying triplet excitations of the two-flavor Schwinger model. For what concerns the 
degenerate triplets of momenta 7r/3 and 5tt/3 (or 2tt/3 and 4tt/3) they do not have definite P-parity, 
but it is always possible to take a linear combination of them with parity opposite to the ground 
state. 



Table 2: Triplet internal quantum numbers 



TRIPLET 


Ai 


A 2 






|1,1,-^^J,0 > 


. /5-2V5 

V 20 


. /5-2V5 

V 20 


1 

2 


1 
2 


|1,1,-S=^.7,0 > 


, / 5+2^ 
V 20 


. / 5+2%/5 
V 20 


3 
2 


3 
2 


|1 1 _5 T 7T > 

1 ) ) 2 ' 3 


8 


8 


3 
2 


1 
2 


|i i 7+VT7 r 2tt ^ 
I 1 ) ± i 4 ' 3 


-2v^-V-2+2VT7 
2+2VT7 


-2 v / 3+V / -2+2VT7 
2+2VT7 


3 
2 


1 

2 


|i i 7+VT7 r 4tt „ 
I 1 -? ± i 4 ' 3 


2v / 3-\/-2+2v'T7 
2+2v / 17 


2v / 3+V- 2 + 2 v'l7 
2+2^17 


1 

2 


3 
2 


|1 1 —5 t 5tt > 


8 


8 


1 

2 


3 
2 



The remaining three triplets in fig. (J^) have no real A's and are characterized by a string of length 
2 and four holes for Q = — |, — |, i, |, i.e. do not belong to the type ( |2.43 ). More precisely, two 
triplets have a string approximately of length 2, due to the finite size of the system, while the triplet 
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of momentum ir has a string exactly of length 2. In table (|^) we summarize the quantum numbers 
of these triplets. 



Table 3: Four holes triplet internal quantum numbers 



TRIPLET 




Ai 


A 2 


IM, 7 -fj, 


2n . 
3 ^ 


2 v / 3-'i\/ 2 + 2 v / T7 
-2+2V17 


2 v / 3+i\/ 2 + 2 v / 17 : 
-2+2\/T7 


|1,1,-J,7T> 




i 
2 


« 
2 




4ir . 
3 ^ 


2V3+iy/ 2+2VI7 


2 v / 3-i-\/ 2 + 2 v / 17 




2-2VT7 


2-2v / 17 



In fig. (j^) it is shown that the spectrum exhibits five singlet states. The lowest lying state 
at momentum 7r is the ground state. Then there are three excited singlets characterized by the 
configuration with two holes ( [2.42] ), i.e. they have one real A and a string of length almost 2. In table 
(^) we summarize their quantum numbers. Among these singlets, those which are not degenerate, 
have P-parity equal to that of the ground state (odd) as it happens in the two-flavor Schwinger 
model. The non- singlet in fact reads 



|0,0,-3J,0 > 



1 



{I mm> +1 titi4t> +1 mm> +i mm> +i utut> +i mm> 
mm> -1 unu> -i mut> -i mm> -i umr>} ■ (2.83) 
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uun> 



The degenerate singlets are again not eigenstates of the P-parity, but it is always possible to take a 
linear combination of them with the a P-parity that coincides with that of the representative state 



(2.83) of the configuration. 



Table 4: Singlet internal quantum numbers 



SINGLET 


A 


As 




|0,0,-3J,0 > 








-1 1 


|0,0,-2J,f > 


V3+2VE 

14 


-2+3V2 

V3(4+v^) 


1 


|0,0,-2J,f > 


\/3+2\/6 


2-3V2 


-1 


14 


V3(4+V2) 



The remaining singlet |0, 0, 



'13 



J, 7r > it is not of the type (2.42). It is characterized by a 



string approximately of length 3 with Ai,i = iy 5+ ?/^ , A24 = and A34 



5+2V13 
12 



Even in finite systems very small like the 4 and 6 sites chains, the "string hypothesis" is a very 
good approximation and it allows us to classify and distinguish among states with the same spin. 

The ground state of the antiferromagnetic Heisenberg chain with N sites is a linear combination 



N 

N 
2 



of all the ( Vr I states with spins up and spins down. These states group themselves into 



14 



E/J* 



tt/3 



2tt/3 



7T 



4tt/3 



5tt/3 



2tt 



-1/2 

-(5- vT3)/2 

-(V- ViV)/2~ 



e 



(5 - v/5)/2 
-3/2 



-5/2 



-e- 



-(7l7 + 7)/4 



-e- 



■o 



-0- 



■G 



-(V5 + 5)/2 



-(v / T3 + 5)/2 



© 

/\ SEVENTHPLET Q TRIPLETS 
<^> QUINTETS P] SINGLETS 



Figure 3: Six sites chain spectrum 



sets with the same coefficient in the linear combination according to the fact that the ground state 
is translationally invariant (with momentum (n) for ^ even (odd)), it is an eigenstate of P-parity 
and it is invariant under the exchange of up with down spins. The states belonging to the same set 
have the same number of domain walls, which ranges from N, for the two Neel states, to 2 for the 
states with ^ adjacent spins up and ^ adjacent spins down. 

The ground state of the 8 sites chain is 



\g.s. >= ^(|V>8 > +a\4 1] > +(3\4 2) > + 7 |4 1) > > > +CI^2 >) (2.84) 

where 



|V>8 >= I 


nnnn>+i mum> 




(2.85) 




> = 


1 TTinni> +1 UUUtT> +translated 


states 


(2.86) 


\4 2) 


> = 


I TTiniTI> +1 IITITTIT> ^translated 


states 


(2.87) 


l# 


> = 


TTIITTI4> ^translated states 




(2.88) 


\4 2) 


> = 


I TTT444T4> +1 II4TTT4T> ^translated 


states 


(2.89) 


\4 3) 


> = 


I TUUm> +1 UTtTUT> ^translated 


states 


(2.90) 


|V>2 >= | 


TTTT4444> ^translated states . 




(2.91) 



By direct diagonalization one gets 

a = -0.412773 (2.92) 
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13 = 0.344301 (2.93) 

7 = 0.226109 (2.94) 

6 = -0.087227 (2.95) 

e = 0.136945 (2.96) 

C = 0.018754 (2.97) 

N = 2 + 16a 2 + 8/3 2 + 4 7 2 + 16<5 2 + 16e 2 + 8C 2 = 6.30356 . (2.98) 

The energy of the ground state is 

Eg. s . = -5.65109J . (2.99) 



Eq.(^99|) differs only by 1.8% from the thermodynamic li mit exp ression E gs . — —8 In 2 = —5.54518 



Moreover also the correlation function of distance 2 Eq. ( [2.101 ) comp uted f or the 8 sites chain is 



G(2) = 0.1957./V, value which is 7% higher than the exact answer Eq.( |2.102D 



In the analysis of finite size systems we were able to find the coefficient (3 of the first set of states 
containing N — 2 domain walls in the ground state. These states are obtained interchanging two 
adjacent spins in the Neel states. The (3 is for a generic chain of JV-sites 

fj= N + 2E g . s . =1 _ 21n2 _ (2100) 
2.3 Spin-spin correlators 

The explicit computation of spin-spin correlations is far from being trivial since the correlator 
G(r) =< g.s.\So ■ S r \g.s. > is not known for arbitrary lattice separations r. For r = 2 it was 
computed by M. Takahashi [^5| in his perturbative analysis of the half filled Hubbard model in 
one dimension. For r > 2 no exact numerical values of G(r) are known. In ]28| ] were given two 
representations of G(r), while in [^9|, [30| the exact asymptotic (r — > oo) expression of G(r) was 
derived. 



In order to explicitly compute the second order energies Eq.( 5.211) and Eq.( 3.212) one has to 
evaluate the correlation function 



1 - 

G(2) = - < 9-s-\S x ■ S x+2 \g.s > (2.101) 

X— 1 



which has been exactly computed in p5| and is given by 

G(2) = -(1 - 161n2 + 9C(3)) = 0.1820 . (2.102) 

In the following we shall show how the knowledge of this correlator allows one to compute explicitly 
the first three "emptiness formation probabilities" , used in Ref. in the study of the Heisenberg 
chain correlators, G{r). The isotropy of the Heisenberg model implies that 

N N 

^2<g.s.\S x -S x+2 \g.s.>=3j2<9^\S 3 x -S 3 x+2 \g.s.> . (2.103) 

x—l x—1 

Let us introduce the probability P3 for finding three adjacent spins in a given position in the 
Heisenberg antiferromagnetic vacuum. Taking advantage of the isotropy of the Heisenberg model 
ground state and of its translational invariance, it is easy to see that the correlator ( 2.103| ) can be 
written in terms of the P^s as 

N 1 

Y, < 9-s.\S 3 x ■ S 3 x+2 \g.s >= N - 2 ( P 3 (TTT) + P 3 (UT) - ^(tU) - ^(4TT) ) ■ (2.104) 

x=l 
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The factor 2 appears in (2.104) due again to the isotropy of the Heisenberg model: the probability 
of a configuration and of the configuration rotated by ir around the chain axis, are the same. 

In pq] the so called "emptiness formation probability" P(x) was introduced. 



P(x)=<g.s.\l[P j \g.s.> , 

3=1 



where 



P; 



1 



(2.105) 



(2.106) 



and Oj is the Pauli matrix. P{x) determines the probability of finding x adjacent spins up in the 
antiferromagnetic vacuum. One gets 

P(ttt) = P(3) 



so that Eq.(2.101) reads 



P(UT) = P{1) ~ 2P(2) + P(3) 
P(ITt) = P(TT|)=P(2)-P(3) 



G(2) = 2P(3) - 2P(2) + -P(l) 



(2.107) 
(2.108) 
(2.109) 

(2.110) 



Using the exact value the correlator G(2) computed in |25| from ( 2.11C ) and from the known values 
of P(2) and P(l) given in §§ 



P(l) 
P(2) 



1 

2 

~(l-ln2) 



(2-111) 

(2.112) 
(2.113) 



one gets 



P(3) = _(l-71n2) + -C(3) 



(2.114) 



For the general emptiness formation probability P(x) of the antiferromagnetic Heisenberg chain, 
an integral representation was given in |2q ], but, to our knowledge, the exact value of P(3) (2.114) 
was not known. 

We illustrate now the computation of a spin-spin correlator which appears in the mass spectrum 
of the two-flavor lattice Schwinger model 

N 

< g.s.\V ■ V\g.s. >= ]T (< g.s.\(S x ■ S y )(S x+1 ■ S y+1 )\g.s. > 

x,y=l 

N 

- < g.s.\(S x ■ S y+ i)(S x+ i ■ S y )\g.s. >) - ^ < g.s.\S x ■ S x+ i\g.s. > . (2.115) 



It is possible to extract a numerical value from Eq.(2.115) only within the random phase approx- 
imation Pq , pi}] . For this purpose it is first convenient to rewrite the unconstrained sum over the 
sites x and y as a sum where all the four spins involved in the VEV's lie on different sites, 



< g.s.\V- V\g.s. > = (< g.s.\(S x ■ S y )(S x+1 ■ S y+1 )\g.s. > 



- < g.s.\(S x ■ S y+1 )(S x+1 ■ S y )\g.s. >) + -N 



\ X! < 9- s -\S x ■ S x+1 \g.s. > < g.s.\S x ■ S x+2 \gs. > (2.116) 



x = l 
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and then factorize the four spin operators in Eq. (2.116) as 

oo 

< g-s.\V ■ V\g.s. > = N^(< g.s.\£> ■ S r \g.s. > 2 - < g.s.\S ■ S r+1 \g.s. >< g.s.\Si ■ S r \g.s. >) 



r=1 



-N ■ 



j N N 

r < 9-s-\S x ■ S x+1 \g.s. > - ^2 < 9- s -\ s x ■ S x+2 \g.s. > (2.117) 



x=l 



Of course, Eq.( 2.117) provides an answer larger than the exact result; terms such as < (...)(...) > 
yield negative contributions which are eliminated once one factorizes them in the form < (...) >< 
(...) >. This is easily checked also by direct computation on finite size systems. 

The spin-spin correlation functions G(r) are exactly known for r — 1,2. For the spin-spin 
correlation functions G(r) up to a distance of r — 30 the results are reported in table (||) |2£| Ej^]. 



Table 5: Spin-spin correlation functions 



r 


G(r) 


r 


G(r) 


1 


-0.4431 


16 


0.0305 


2 


0.1821 


17 


-0.0296 


3 


-0.1510 


18 


0.0274 


4 


0.1038 


19 


-0.0267 


5 


-0.0925 


20 


0.0249 


6 


0.0731 


21 


-0.0242 


7 


-0.0671 


22 


0.0228 


8 


0.0567 


23 


-0.0223 


9 


-0.0532 


24 


0.0211 


10 


0.0465 


25 


-0.0206 


11 


-0.0442 


26 


0.0196 


12 


0.0395 


27 


-0.0193 


13 


-0.0379 


28 


0.0183 


14 


0.0344 


29 


-0.0181 


15 


-0.0332 


30 


0.0172 



For r > 30, one may write |2J 



G(r) 



with g(r) satisfying 



and 



- - 1 n 3 

4 V^ 3 r^W 16 
^((-l)- + l + 0(ff(r)) + 0(i)] 



7(r) - ^r{l + \g(r)Hg(r))) 



C(r) = ln(2V2^e 7+1 r) 



Eq.(2.119) may be solved by iteration. To the lowest order in ^ one finds 

1 1 

g(r) 



C(r) C(r) : 



TnC(r) 



(2.118) 

(2.119) 
(2.120) 

(2.121) 
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Inserting fl2.12l|) in Eq.(|2.118|) leads to 



G {r)~^\VW)[l + ^^C{r)) + 0{-^) . (2.122) 



Inserting Eq. fl2.122 ) in (2.117), one finally gets 

< g.s.\V ■ V\g.s. >= 0.38167V (2.123) 

Last, we report the following exact three spin correlators that have been used in the determi- 
nation of the mass spectrum of the two-flavor lattice Schwinger model 

N 



<g.s.\J2S^ +1 S 3 x+2 \g.s.> = (2.124) 

x=l 
N 

<g.s.\J2siS- +1 S 3 x+2 \g.s.> = (2.125) 



x = l 

N 



<g.s.\J2S-S+ +1 S x+2 \g.s.> = (2.126) 

x=l 
N 

<g.s.\J2SiS 3 +1 S- +2 \g.s.> = (2.127) 



x = l 
N 



<g.s.\J2S-S x+1 S+ +2 \g.s.> = (2.128) 



x=\ 



So, on a spin singlet, not only the VEV of J2x=i * s zero i but also every VEV with an odd 
number of S 3 . 

2.4 SU(Af) quantum antiferromagnetic chains 

It is our purpose to introduce now spin- 1/2 antiferromagnetic Hcisenberg chains where "spins" are 
generators of the SU(Af) group. In the limit J\f = 2 one has the usual antiferromagnetic Heisenberg 
chain discussed in the previous section. An U(Af) spin-1/2 quantum antiferromagnetic chain is 
described by the Hamiltonian 

N 

H U j W = jY J S ab {x)S ba {x + l) (2.129) 

x=l 

where S ab (x) a,b= 1, . . . , N are the generators of U(N) satisfying the Lie algebra 

[S ab (x), S cd {y)} = {S ad {x)S bc - S cb (x)8 ad )8 xy (2.130) 
and they can be conventionally represented by fermion bilinear operators 

Sab {x)=^ bx - 5 -f . (2.131) 

The representation of the algebra on each site is fixed by specifying the fermion number occupation 

M 

P {x)=Y,Saa{x) . (2.132) 

By fulfilling the global neutrality condition Ylx=i P( x ) = 0> one ma y choose 

M r 

m x even 



a=l 
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or viceversa the opposite choice for x even-odd. Eq.(2.133) restricts on each site the Hilbert space 
to a representation with Young tableau of m rows for x even and Af — m rows for x odd. For each 
site x p(x) is the generator of the U(l) subgroup of U(Af). 

Let us use the basis T a = (T Q )*, a = 1, . . . ,7V 2 - 1, of the Lie algebra of SU(Af) in the 
fundamental representation such that tr(T a T^) = 6 a/3 /2 and [T a , T^] = if^T*, where f a ^ are 
the structure constants. By means of 



1 ab 1 cd 



■^SadSbc — -^Jj5ab8cd 



and redefining the group generators 



s° x = to 



one can rewrite the Hamiltonian (2.129) as 

AT 



1) + H 



SU(Af) 
J 



(2.134) 



(2.135) 



(2.136) 



where 



H 



N 

SU(Af) _ j ST^ qa qo 

j — J 2.^1 x x 

x=l 



(2.137) 



is the Hamiltonian of an SU(Af) quantum antiferromagnet. From Eq.( ^.136 ) it is clear that by 
fixing p(x) is reduced to Hj U ^\ 

The representations of the SU(Af) spins relevant for the relationship of the spin Hamiltonians 



( 2.137 ) with the Schwinger models are different when Af is even or odd. When Af is even we shall 
consider the representation on each site with Young tableau of one column and Af/2 rows so that 
p(x) = 0. When Af is odd we shall consider on one sublattice the representation with Young tableau 
of one column and (Af + l)/2 rows and of one column and (Af — l)/2 rows on the other sublattice. 

As an exp licit example let us consider the case Af = 3. The generators T a of the SU(3) group 
in Eq.( ^.135 ) are the Gell Mann matrices. By denoting u,d and s the three flavors, the eight spin 
operators S" read 



&x 




(2.138) 


S. j. 


= -i^lx^dx + W d J>ux 


(2.139) 




= i>L*Pux - ^teipdx 


(2.140) 


S x 


= ^ixi>sx + lft x lfux 


(2.141) 


s 5 x 


= -i^lx^sx + ilptx^ux 


(2.142) 


&x 


= l/jfolpsx + IpLlpdx 


(2.143) 




= -Wdx^™ + ^Ix^dx 


(2.144) 




= -j^lx^ux + i>dJ>dx - 2V>L^x) 


(2.145) 



Let us find the ground state of the Hamiltonian (2.137) for a chain of two sites. Taking the 
representations with one particle on site 1 and two particles on site 2, the ground state with energy 



E, 



g.s. 



-16 J/3 reads 



\g.s. >= -j=(\u d s > -\d 



> 



>) 



(2.146) 



The state ( [2.146 ) is a singlet of SU(3), i.e. it is annihilated by the Casimir S 2 = (Si + S2) 2 = 
Sa=i [(^l*) 2 + (S2) 2 + 25f S<f\ . If one chooses the repre sentat ion with two particles on site 1 and 
one particle on site 2, the ground state degenerate with ( [2. 146 ) is 



>- T3 ( 



I Lb , I Lb x 

u> -| d>+\ , s>) 



(2.147) 
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Diagonalizing the translation operator T with the Hamiltonian ( [2.137 ), one gets two degenerate 
ground states 

>±= \g.s. > ±\g.s. >' 

V2 



In Eq.( p.l48 ) the linear combination with the + has momentum zero and the one with the — has 



momentum tt. By studying this very simple example, one can infer that the thermodynamic limit 
N — > oo analysis for a generic SU(Af) chain is very involved. 

Unfortunately, no analysis with a level of completeness such as the one given in |18j for the 
SU(2) case does exist for a generic symmetry group SU(Af) with spins in the representation in 
which we are interested. In [p3| SU (Af) antiferromagnetic models were solved for a particular spin 



representation such that the Hamiltonian ( p.l37| ) becomes 



N 

rSC/(AO 



H suW =J J2P xx+1 (2.149) 

x=l 

with P X x+i the operator that permutes whatever objects occupy sites x and x + 1. The spectrum 



of Eq.(2.149) was shown to exhibit massless excitations. 

A large Af expansion approach has been performed in j3^] for an SU (Af) antiferromagnetic chain 
characterized by spins living in a representation with Young tableaux of one row on one sublattice 
and Af — 1 rows on the other sublattice. In the case of Young tableaux of one column it was 
found that the ground state is twofold degenerate and breaks translational and parity symmetry. 
Moreover the elementary excitations are massive non relativistic solitons in the large Af limit with 
a mass of 0(Af). 

In J35| the Lieb-Shultz-Mattis theorem |36) was generalized to SU(Af) spin chains. The theorem 
proves that a half-integer-S spin chain with essentially any reasonably local Hamiltonian respecting 
translational and rotational symmetry either has zero gap or else has degenerate ground states 
spontaneously breaking translational and parity invariance. In |36| ] it was proved the existence of a 
unique ground state of the SU(Af) chains where Af is even for spins in the antisymmetric Af tensor 
representation. Under the assumption of a unique ground state which must be a SU(Af) singlet, 
an infinitesimal energy gap was found for all the representations of SU (Af) whose Young tableaux 
contain a number of boxes not divisible by Af. Of course the degeneracy of the ground state trivially 
implyes a zero gap. 

Using the Lieb-Shultz-Mattis theorem we can state that when Af is even the ground state of the 



Hamiltonian (2.137) is unique and there are gapless excitations for the spin representation with a 
Young tableau of one column and Af/2 rows. When Af is odd the ground state is twofold degenerate, 
as we illustrated with the simple example of the two site SU(3) chain. It is doubtful if there exist 
gapless excitations in this case and at present time we are investigating this problem. 



3 The two-flavor lattice Schwinger model 

The one-flavor lattice Schwinger model has been studied in |t4| . The solution of the strong coupling 
problem is identical to solving a particular type of Ising system with long range interaction. We 
showed that the mass of the elementary excitations and the chiral condensate could be computed 
reliably from an extrapolation to weak coupling using Pade approximants. Our analysis improved 
previous ones ]37j by taking careful account of all discrete symmetries of the continuum theory. We 
do not discuss here the one-flavor model, since it is mapped into the classical Ising spin model and 
it is out of the scope of these lectures. 

We study the ST/(2)-fiavor lattice Schwinger model in the hamiltonian formalism using staggered 
fermions. The existence of the continuum internal isospin symmetry makes the model much more 
interesting than the one-flavor case; the spectrum is extremely richer, exhibiting also massless 
excitations and the chiral symmetry breaking pattern is completely different from the one-flavor 
case. We shall demonstrate 113, 03] that the strong coupling limit of the two-flavor lattice Schwinger 
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model is mapped onto an interesting quantum spin model — the one-dimensional spin-1/2 quantum 
Heisenberg antiferromagnet. The ground state of the antifcrromagnetic chain has been known since 
many years [ p~7| and its energy was computed in p3fl ; the complete spectrum as been determined 
by Faddeev and Takhtadzhyan Q using the algebraic Bethe ansatz. 

The two-flavor lattice Schwinger model with non-zero fermion mass m has been analysed in ]38| ] 
in the limit of heavy fermions m 3> e 2 ; good agreement with the continuum theory has been found. 

There are by now many hints at a correspondence between quantized gauge theories and quantum 
spin models, aimed at analyzing new phases relevant for condensed matter systems |^|, ^, f^, |8|, |9| . 
Recently Laughlin has argued that there is an analogy between the spectral data of gauge theories 
and strongly correlated electron systems Moreover, certain spin ladders have been shown to 
be related to the two-flavor Schwinger model Q. 

The correspondence between the SU(2) flavor Schwinger model and the quantum Heisenberg 
antifcrromagnetic chain provides a concrete computational scheme in which the issue of the cor- 
respondence between quantized gauge theories and quantum spin models may be investigated. 
Because of dimensionality of the coupling constant in (l-fl)-dimcnsions the infrared behavior is 
governed by the strong coupling limit, and it is tempting to conjecture the existence of an exact 
correspondence between the infrared limits of the Heisenberg and two-flavor Schwinger models. 
We shall derive |l^, |l5| results which support this conjecture. For example the gapless modes in 
the spectra have identical quantum numbers; within the accurancy of the strong coupling limit, 
the gapped mode of the two-flavor Schwinger model was also identified in the spectrum of the 
Heisenberg model. 

In this section we present a complete study Jl9|, [D| |l6) of the strong coupling limit of the 
two-flavor lattice Schwinger model. We firstly compute explicitly the masses of the excitations to 
the second order in the strong coupling expansion; this computation needs the knowledge of the 
spin-spin correlators of the quantum Heisenberg antiferromagnetic chain. The continuum massless 
two-flavor Schwinger model does exhibit neither an isoscalar (ipip) nor an isovector ^cr ?/;^ chiral 
condensate, since this is forbidden by the Coleman theorem p2fl . On the lattice these fermion 
condensates are zero to all the orders in the strong coupling expansion pl| . The pertinent non-zero 

chiral condensate is < ij^L ^^r4^r > an d we computed its lattice expression up to the second 
order in the strong coupling expansion Q. It should be noticed that, in abscence of gauge fields, 
the chiral condensate is zero, is different from zero only when the fermions are coupled to gauge 
fields. This can be viewed as the manifestation of the chiral anomaly in this model. 



3.1 The model 



The action of the 1 + 1-dimensional electrodynamics with two charged Dirac spinor fields is 

2 



S = jd 2 x ^^j^ + ^A^a-^F^ 
\_a=l c 



(3.150) 



The theory has an internal SUl(2) (g> S , t r ^(2)-flavor isospin symmetry; the Dirac fields are an 
isodoublet whereas the electromagnetic field is an isosinglet. It is well known that in 1 + 1 dimensions 
there is no spontaneous breakdown of continuous internal symmetries, unless there are anomalies 
or the Higgs phenomenon occurs. Neither mechanism is possible in the two- flavor Schwinger model 
for the SUl(2) <g> 5t/_R,(2)-symmetry: isovector currents do not develop anomalies and there are no 
gauge fields coupled to the isospin currents. The particles belong then to isospin multiplets. For 
what concerns the U(l) gauge symmetry there is an Higgs phenomenon |il|j . 

The action is invariant under the symmetry 



SU L {2) ® SU R (2) <g> U v (l) ® U A {1) 
The group generators act on the fermion isodoublet to give 



(3.151) 



SU L {2) : ip a (x) — > (e^ — ^) ab i> b {x) , x/, a (x) — » ^ h {x)(e 



<^Pr 



(3.152) 
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SU R (2) 






) , i>a(x) 




Mi) 




~> (e^U ^(x) 




-VjC*) (e" ie(a:31 )6a 


Mi) 


4>a(x) - 




lib) - 





where cr Q are the Pauli matrices, 9 a , 9{x) and a are real coefficients and 



P L = - 2 (l 



75 j 



2(1+75) 



(3.153) 
(3.154) 
(3.155) 



(3.156) 



At the classical level the symmetries (3.152— ft. 155 ) lead to conservation laws for the isovector, 
vector and axial currents 



tW7^(y)«A(^) 

^a( x h >11 ab1pb(x) 



(3.157) 
(3.158) 

(3.159) 
(3.160) 



It is well known that at the quantum level the vector and axial currents cannot be simultaneously 
conserved, due to the anomaly phenomenon [ |43"| . If the regularization is gauge invariant, so that 
the vector current is conserved, then the axial current acquires the anomaly which breaks the 
V A ( 1 )-symmetry 



WW = 2^F^(x) 



(3.161) 



The isoscalar and isovector chiral condensates are zero due to the Coleman theorem p^] ; in fact, 
they would break not only the Ua(X) symmetry of the action, but also the continuum internal 
symmetry SU L {2) <g> SU R {2) down to SU V (2). There is, however, a SU L {2) ® SU R (2) invariant 
operator, which is non-invariant under the J7 J 4(l)-symmetry; it can acquire a non- vanishing VEV 
without violating Coleman's theorem and consequently may be regarded as a good order parameter 
for the J7 J 4(l)-breaking. Its expectation value is given by 



42 



47T 7T 



(3.162) 



It describes a process in which two right movers are anihilated and two left movers are created. 
Note that F, being quadrilincar in the fields, is actually invariant under chiral rotations of ir/2, 
namely under the discrete axial symmetry 



4>a(x) -> ~f 5 4>a(x) 1p a {x) 



(3.163) 



As a consequence, this part of the chiral symmetry group is not broken by the non-vanishing VEV 
of F ( p62| ). 

The lattice theory faithfully reproduces the pattern of symmetry breaking of the continuum 
theory; this happens even if on the lattice the S'f7(2)-flavor symmetry is not protected by the 
Coleman theorem. The isoscalar and isovector chiral condensates are zero also on the lattice, 
whereas the operator F acquires a non- vanishing VEV due to the coupling of left and right movers 
induced by the gauge field. The continuo us axi al symmetry is broken explicitly by the staggered 
fermion, but the discrete axial symmetry (3.163) remains. 

The action ( 3.15C ) may be presented in the usual abelian bosonized form [ft7| . Setting 

: ^ al ^ a ■= -U^d^ , a = 1, 2 , (3.164) 
the electric charge density and the action read 

jo =: V4V>1 + V4V>2 := + $ 2 ) (3.165) 
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s 



d 2 x 



Z I Z7T 



By changing the variables to 



(3.166) 



-^($l+$ 2 ) 



$_ = ($ t _$ 2 ) 



one has 



5 = / d 2 x f + -du$-d»$- - ^<f> 2 

J \2 2 7T 



(3.167) 
(3.168) 

(3.169) 



The theory describes two scalar fields, one massive and one massless. $ + is an isosinglet as 



evidenced from Eq.(3.165); its mass ms 



e c comes from the anomaly Eq.(3.161) ED]. Local 



electric charge conservation is spontaneously broken, but no Goldstone boson appears because the 
Goldstone mode may be gauged away. represents an isotriplet; it has rather involved nonlinear 
transformation properties under a general isospin transformation. All three isospin currents can be 
written in terms of but only the third component has a simple representation in terms of $_; 
namely 



(3.170) 



The other two isospin currents j^{x) and j 2 (x) are nonlinear and nonlocal functions of $_ p7| ; a 
more symmetrical treatment of the bosonized form of the isotriplet currents is available within the 
framework of non abelian bosonization p8[ | . For the multiflavor Schwinger model this approach has 
been carried out in , providing results in agreement with jlTj . 

The excitations are most conveniently classified in terms of the quantum numbers of P-parity 
and G-parity; G-parity is related to the charge conjugation G by 



G 



$_ is a G-even pseudoscalar, while $+ is a G-odd pseudoscalar 



I PG = 1- 

jPG = Q - 



(3.171) 



(3.172) 
(3.173) 



The massive meson $ + is stable by G conservation since the action (3.169) is invariant under 
$ + — ► -$+. 

In the massive SU (2) Schwinger model — when the mass of the fermion m is small compared to 
e 2 (strong coupling) — Coleman j|7j showed that - in addition to the triplet <f>_ (I PG = 1 h ) the 
low-energy spectrum exhibits a singlet I PG = ++ lying on top of the triplet <&_. In this limit the 
gauge theory is mapped to a sine-Gordon model and the low lying excitations are soliton-antisoliton 
states. When m — > 0, these soliton-antisoliton states become massless [pp| ; in this limit, the analysis 
of the many body wave functions, carried out in ref. |Eo|, hints to the existence of a whole class of 
massless states with positive G-parity; P-parity however cannot be determined with the procedure 
developed in These are not the only excitations of the model: way up in mass there is the 

isosinglet I PG = , (<&+), already discussed in ref. |l7|. The model exhibits also triplets, whose 
mass — of order mj or greater — stays finite p3] ; among the triplets there is a G-even state []. 

2 K. Harada private communication. 
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The Hamiltonian, gauge constraint and non-vanishing (anti-) commutators of the continuum 
two-flavor Schwinger model are 



H 



dx 



^E\x) + ELi i>Ux)a {id x + eA(x)) $ a (xj 

d x E( X ) +Ea=lV>i(2#a(*)~0 

[A{x),E(y)] = i6{x-y) , [^ a {x), ^ (*/)} = S ab S(x - y) 



(3.174) 

(3.175) 
(3.176) 



A lattice Hamiltonian, constraint and (anti-) commutators reducing to (3.174,3.175, ^.176 ) in the 
naive continuum limit are 



, N 

e a 



x — 1 x—1 a—1 

E x -E x -i + iA^hx + ^2,x^2.x ~ 1 ~ , 
[A x ,E y ] = i6 x , y , |^ a ,x, M>,y\ = $abSxy ■ 



N 2 



x + l 



(3.177) 



The fermion fields are defined on the sites, x = 1, ...,N, the gauge and electric fields, A x and E x , 
on the links [a;; x+l], N is an even integer and, when TV is finite it is convenient to impose periodic 
boundary conditions. When N is finite, the continuu m lim it is the two-flavor Schwinger model on a 
circle [Q. The coefficient t of the hopping term in ( [3.177 ) plays the role of the lattice light speed. 
In the naive continuum limit, = e c and t = 1. 

The Hamiltonian and gauge constraint exhibit the discrete symmetries 
• Parity P: 

A x — > —A-x-i, E x — > i> a , x — (-1)^,-*, V>L — » (-IJ^Lx (3-178) 



• Discrete axial symmetry T: 

► A x+ i, E x — 

• Charge conjugation C: 

A x ► Ax+i) E x > - 

• G-parity: 

A x 

4>i,x 

Tp2,x ► - 



E x+ 1, 1p a , x ► lpa.x+1, 1pl,x > V4,x+1 (3.179) 

E x+1 , lp a . x ► % x+1 , i>l x ► tl>a,x+l (3.180) 



— ► — A x +i, E x — > — E x+ i 

— 4»+i> — » (3.181) 



The lattice two-flavor Schwinger model is equivalent to a one dimensional quantum Coulomb 
gas on the lattice with two kinds of particles. To see this one can fix the gauge, A x = A (Coulomb 
gauge). Eliminating the non-constant electric field and using the gauge constraint, one obtains the 
effective Hamiltonian 



H s = H u + H p 
it 



x,y 



+ 



2a 



J2&l*+i eiA ^>* - <* e ~ iA v>a,x+o 



where the charge density is 



P( x ) = V'UV'l.a: + ^tx^x - 1 , 



(3.182) 
(3.183) 
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and the potential 

1 N ^ 1 

V(x - y) = — V e «*n(«.-»)/W ( 3 . 184 ) 

v * ; iV ^ 4sin 2 ^ V y 

is the Fourier transform of the inverse laplacian on the lattice for non zero momentum. The constant 
modes of the gauge field decouple in the thermodynamic limit N — ► oo. 



3.2 The strong coupling limit and the antiferromagnetic Heisenberg Hamil- 
tonian 



In the thermodynamic limit the Schwinger Hamiltonian (3.182), rescaled by the factor e^a/2, reads 

H = H + eH h (3.185) 

with 

'{x-yf 



Ho 



E 



TV 

x>y 

—i(R — L) 



- y) 



p{x)p{y) 



(3.186) 
(3.187) 



and e = t/e 2 L a 2 . In Eq.(3.187) the right R and left L hopping operators are defined (L = RJ) as 



N 



N 



N 



(3.188) 



x—l x—1 a—1 x—1 a— 1 

On a periodic chain the commutation relation 

[R,L]=0 (3.189) 

is satisfied. 

We shall consider the strong coupling perturbative expansion where the Coulomb Hamiltonian 
(3.186) is the unperturbed Hamiltonian and the hopping Hamiltonian (3.187) the perturbation. 
Due to Eq.(3.183) every configuration with one particle per site has zero energy, so that the ground 
state of the Coulomb Hamiltonian ( 3.186| ) is 2^ times degenerate. The degeneracy of the ground 
state can be removed only at the second perturbative order since the first order is trivially zero. 

At the second order the lattice gauge theory is effectively described by the antiferromagnetic 
Heisenberg Hamiltonian. The vacuum energy — at order e 2 — reads 



E, 



(2) 



=< H 



n 



&a — no 



H h > 



(3.190) 



where the expectation values are defined on the degenerate subspace and n is the operator projecting 
on a set orthogonal to the states with one particle per site. Due to the vanishing of the charge 
density on the ground states of i?o, the commutator 



[Ho, Hh] = Hh 



(3.191) 



holds on any linear combination of the degenerate ground states. Consequently, from Eq.( [3.19c| ) 
one finds 



E, 



(2) 



-2 < RL > 



(3.192) 



On the ground state the combination RL can be written in terms of the Heisenberg Hamiltonian. 
By introducing the Schwinger spin operators 



Sx = ipt, x -^-ipb, t 



(3.193) 
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the Heisenberg Hamiltonian Hj reads 



N 



x=l 

and, on the degenerate subspace, one has 

N 



Hj — ( S x ■ S x+ i — 

x = l ^ 

N ( 1 I 

= Y.\-2 LxRx ~\ p[x)p{x + X) ) (3 - 194) 



<Hj>=(j2(s x -S x+1 -\) ) = {J2(-Il x R x ) ) . (3.195) 



\x=l 



Taking into account that products of L x and R v at differe nt points have vanishing expectation 
values on the ground states, and using Eq.( 3.195 ), Eq.( 3.192] ) reads 



E<g ] = 4<Hj> . (3.196) 

The ground state of Hj singles out the correct vacuum, on which to perform the perturbative 
expansion. In one dimension Hj is exactly diagonalizable [ p7| 20 . In the spin model a flavor 1 



particle on a site can be represented by a spin up, a flavor 2 particle by a spin down. The spectrum of 
Hj exhibits 2 N eigenstates; among these, the spin singlet with lowest energy is the non degenerate 
ground state \g.s. >. 

We shall construct the strong coupling perturbation theory of the two-flavor Schwinger model 
using \g.s. > as the unperturbed ground state. \g.s. > is invariant under translations by one lattice 
site, which amounts to invariance under discrete chiral transformations. As a consequence, at 
variance with the one-flavor model ]14| ], chiral symmetry cannot be spontaneously broken even in 
the infinite coupling limit. 

\g.s. > has zero charge density on each site and zero electric flux on each link 

p{x)\g.s. >= , E x \g.s. >= (x=l,...,N) . (3.197) 

\g.s. > is a linear combination of all the possible states with ^ spins up and ^ spins down. The 
coefficients are not explicitly known for general N. In section 2.2 we exhibited \g.s. > explicitly for 
finite size systems of 4, 6 and 8 sites. The Heisenberg energy of \g.s. > is known exactly and, in 
the thermodynamic limit, is |2^, [ll| 

Hj\g.s. >= (-JV In 2)\g.s. > . (3.198) 
Eq.( 3.198 ) provides the second order correction Eq. ( 3.196] ) to the vacuum energy, Efl = -AN In 2. 



There exist two kinds of excitations created from \g.s. >; one kind involves only spin flipping 
and has lower energy since no electric flux is created, the other involves fermion transport besides 
spin flipping and thus has a higher energy. For the latter excitations the energy is proportional to 
the coupling times the length of the electric flux: the lowest energy is achieved when the fermion is 
transported by one lattice spacing. Of course only the excitations of the first kind can be mapped 
into states of the Heisenberg model. 

In Jl8| the antiferromagnetic Heisenberg model excitations have been classified. There it was 
shown that any excitation may be regarded as the scattering state of quasiparticles of spin-1/2: 
every physical state contains an even number of quasiparticles and the spectrum exhibits only 
integer spin states. The two simplest excitations of lowest energy in the thermodynamic limit are 
a triplet and a singlet Jl8| ; they have a dispersion relation depending on the momenta of the two 
quasiparticles. For vanishing total momentum (relative to the ground state momentum P g . s . = 
for y even, P g . s . — n for -y odd) in the thermodynamic limit they are degenerate with the ground 
state. 

In section 2.2 we showed that even for finite size systems, the excited states can be grouped 
in families corresponding to the classification given in jl8|. We explicitly exhibited all the energy 
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eigenstates for TV — 4 and N = 6. The lowest lying are a triplet and a singlet, respectively; they 
have a well defined relative (to the ground state) P-parity and G-parity - 1 h for the triplet and 
++ for the singlet. Since they share the same quantum numbers these states can be identified, in 
the limit of vanishing fermion mass, with the soliton-antisoliton excitations found by Coleman in 
his analysis of the two-flavor Schwinger model. A related analysis about the parity of the lowest 
lying states in finite size Heisenberg chains, has been given in 

Moreover in a whole class — A4 af — of gapless excitations at zero momentum was singled out 
in the thermodynamic limit; these states are eigenstates of the total momentum and consequently 
have positive G-parity at zero momentum. The low lying states of the Schwinger model also contain 
|fio| many massless excitations with positive G-parity; they are identified ^] with the excitations 
belonging to M.af- The mass of these states in the Schwinger model can be obtained from the 
differences between the excitation energies at zero momentum and the ground state energy. The 
energies of the states \ex. > belonging to the class M.af have the same perturbative expansion of 
the ground state. Consequently, the states \ex. > at zero momentum up to the second order in the 
strong coupling expansion have the same energy of the ground state ( |3.190| ), Eex — —47V In 2. To 
this order the mass gap is zero. Higher order corrections may give a mass gap. 



3.3 The meson masses 



In this section we determine the masses for the states obtained by fermion transport of one site 
on the Heisenberg model ground state. Our analysis shows that besides the G-odd pseudoscalar 
isosinglet with mass ms — eLy/Z/n, there are also a G-even scalar isosinglet ++ and a 
pseudoscalar isotriplet 1 h and a G-odd scalar isotriplet 1 H with masses of the order of ms or 
greater. The quantum numbers are relative to those of the ground state Ig G = ++ for N/2 even 
Iq.?. = 0~~ for N/2 odd. 



Two states can be created using the spatial component of the vector j 1 (x) Eq.(3.159) and 
isovector j^(x) Eqs.(3.157, 3. 15S ) Schwinger model currents. They are the G-odd pseudoscalar 
isosinglet I PG = and the G-even pseudoscalar isotriplet I PG = 1 h . The lattice operators 
with the correct quantum numbers creating these states at zero momentum, when acting on \g.s. >, 
read 



A? 



S = R + L = J2jHx) 



(3.199) 



TV 



T+ = (T_)t = i?< 12 ) + L< 12 ) =^j|(a 



x=l 



N 



T = -^(i? (11) +l( 11 >~rW -L^)^J2jI(x) 
V 2 



(3.200) 



(3.201) 



i?( afc ) and 7> fc ) in ( ^2^ , ^20l| ) are the right and left flavor changing hopping operators (£A ' 
(i?( afc ))t) 



JV 



The states are given by 



They are normalized as 

< S\S > 
<T+\T+ > 



\S> = \Q—>=S\g.s.> (3.202) 

\T ± > = |1 _+ ,±1 >=T±\g.s. > (3.203) 

|T > = |1-+ O>=T o |0.s. > . (3.204) 

< g.s.\S^S\g.s. >= -4 < g.s.\H,,\g.s. >= 4jVln2 (3.205) 

^(N+ < g.s.\Hj\g.s. >) = ^jV(1 - In 2) (3.206) 
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and 



< T \T >=< T_|T_ >=< T+\T+ > 



(3.207) 



In Eqs. (3. 205, 3. 206, 3.207) < g.s.\g.s. >= 1. 



The isosinglet energy, up to the second order in the strong coupling expansion, is Es = E s ^ 



e 2 E i g ) with 



E 



(0) < S\H \S > 

< S\S > 
: Ji < S\HjA s H h \S > 



Ac 



JLs_ 



E » <s\s> 

and 1 — lis a projection operator onto \S >. On \g.s. > 
[H , (U s H h ) n S] = (n + l)(U s H h ) n , (n = 0, 1, . 



holds; Eq. (3.209) may then be written in terms of spin correlators as 



Ex=l < 9-S-\S x ■ S : 



E 



(2) 



x+2 



\g.s. > 



< g-s.\Hj\g.s. > 



(3.208) 
(3.209) 

(3.210) 
(3.211) 



One immediately recognizes that the excitation spectrum is determined once < g.s.\S x ■ S x +2\g-S- > 
is known. Equations similar to Eq. ( ft. 211 ) may be established also at a generic order of the strong 
coupling expansion. 

At the zeroth perturbative order the pseudoscalar triplet is degenerate with the isosing let EP = 



Eg = 1. Following the same procedure as before one may compute the energy of the states (3.203) 



and ( 3.204 ) to the second order in the strong coupling expansion. To this order, the energy is given 

by 



E (2) - E {2) 



A DS (T) 



4 < g.s.\Hj\g.s. > -r,V v , < g . s .\S x ■ S : 



x+2 



\g.s. > 



N+ < g.s.\Hj\g.s. > 

where in terms of the vector operator V — Ylx=i ^ x ^ ^+1, one can write Ads(T) as 

,< g.s.KVrflg.s. > + < g.s.\{V 2 ) 2 \g.s. > 



Ads(T±) 
Ads (T ) 



= 12 



= 12 



N+ < g.s.\Hj\g.s. > 
2<g.s.\(V 3 ) 2 \g.s.> 



N+ < g.s.\Hj\g.s. > 



(3.212) 

(3.213) 
(3.214) 



The VEV of each squared component of V on the rotationally invariant singlet \g.s. > give the 
same contribution i.e. Aos(T±) = Ads(Tq): the triplet states (as in the continuum theory) have a 
degenerate mass gap. This is easily verified by direct computation on finite size systems; when the 
size of the system is finite one may also show that A us is of zeroth order in N. 



The excitation masses are given by ms = ^rr(Es — E g , s .) and ttit 



''I." 

2 



(E T — Eg. s ). Conse- 



(2") (2) 

qucntly, the (A^-dependent) ground state energy terms appearing in Eg and E T cancel and what 
is left are only N independent terms. This is a good check of our computation, being the mass an 
intensive quantity. 

In principle one should expect also excitations created acting on \g.s. > with the chiral currents, 
in analogy with the one flavor Schwinger model where, as shown in ref. |j3~7f , the chiral current 
creates a two-meson bound state. The chiral currents operators for the two flavor Schwinger model 
are given by 



(3.215) 
(3.216) 
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The corresponding lattice operators at zero momentum are 

s 5 =#-£ = E"=ii 5 (*) 



= (T 5 )t=jR (12)_ jL (12) =E ^ ij 5 (a;) 

The states created by ( 3.217|]3~218| , |3.219 ) when acting on \g.s. >, are 

\S 5 > = |0++ >= S 5 |g.s. > 
|T|> = |l+-,±l>=T|| 5 .s. > 
|T 5 > = |1+-,0>=T 5 | 5 . S . > . 



They are normalized as 

< S 5 \S 5 > 

< Tl\Tl > 

and 



< g.s.\S 5 ^S 5 \g.s. >= -4 < g.s.\Hj\g.s. >= 4iVlog2 
^(7V+ < g.s.\Hj\g.s. >) = |iV(l - log 2) 



< T^T* >=< T^T 5 >=< T||T| > 



(3.217) 
(3.218) 
(3.219) 



(3.220) 
(3.221) 
(3.222) 



(3.223) 
(3.224) 

(3.225) 



Following the computational scheme used to study \S > and \T >, one finds for the state \S 5 > 



E 



E 



(0) 

S 5 

(2) 
S 5 



1 



E <2) 

as. 



12-3 



£s=i < 9-s-\S x ■ S x+2 - \ \g.s. > 
< 9-S-\Hj\g.s. > 



For the triplet |T 5 > one gets 



E, 



(0) 



(2) _ 

^5 



1 



g-s- ~ 



E I= i < 9-s.\S x ■ S x+2 - \\g.s. > -4 < g.s.\Hj\g.s. > 



(3.226) 
(3.227) 

(3.228) 
(3.229) 



N+ < g.s.\Hj\g.s. > 

Now we can compute the mass spectrum up to the second order in the strong coupling expansion 



Using Eq.(2.102), the isosinglet mass reads as 

1 



^ = - + 1.9509 e 2 

e A a 2 



For what concerns the isotriplet mass, since the double sum in Eq.( 3.212] ) is given by 

< g.s.\V ■ V \g.s. > 



(T) 



N+ < g.s.\Hj\g.s. > 



using Eq.(2.123), one gets 



1^ = 1 +0.0972 e 2 

e z L a 2 



(3.230) 



(3.231) 



(3.232) 



The existence of massive isotriplets was already noticed in ]50| , and their mass in the continuum 
theory was numerically computed for various values of the fermion mass. In particular there is a 
G-parity even isotriplet with mass approximately equal to the mass of the isosinglet . 

The masses of the jS's > isosinglet and the \T§ > isotriplet are 



m S 5 



m T 5 



5.85e 2 
4.4069e 2 



(3.233) 
(3.234) 
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Equations ( |3.230 ), (3.232), ( 3.233| ) and (3.234) provide the values of ms, mr, ms 5 and rn T r, for 
small values of z = e 2 — -4—i up to the second order in the strong coupling expansion. Whereas 

( 3.232 ) is only approximate ( 3.230 ), ( 3.233 ) and ( 3.234 ) are exact at the second order in the e 
expansion. We extrapolated these masses to the continuum limit using the standard technique of 
the Pade approximants |iq] and we got results in good agreement with the continuum. 



4 The multiflavor lattice Schwinger models 

We now study the TV-flavor lattice Schwinger models in the hamiltonian formalism using staggered 
fermions. The TV-flavor Schwinger models have many features in in common with four dimensional 
QCD: at the classical level they have a symmetry group U L {M) ® U R (N) = SU L (N) ® SU R (J\f) (g) 
Uv(i) (8 UaW that is broken down to SUl{N) ® SUr{N) ® tV(l) by the axial anomaly exactly 
like in QCD [^2[. The massless TV-flavor Schwinger models describe no real interactions between 
their particles as one can infer by writing the model action in a bosonizcd form. The model exhibits 
one massive and TV 2 — 1 massless pseudoscalar "mesons" J53| . 

On the lattice we shall prove that — at the second order in the strong coupling expansion — 
the lattice Schwinger models are effectively described by SU(N) quantum antiferromagnetic spin- 
1/2 Heisenbcrg Hamiltonians with spins in a particular fundamental representation of the SU(Af) 
Lie algebra and with nearest neighbours couplings. The features of the model are very different 
depending on if TV is odd or even. When TV is odd, the ground state energy in the strong coupling 
limit is proportional to e\, the square of the electromagnetic coupling constant. In contrast, when 
TV is even the ground state energy in the strong coupling limit is of order 1. This difference arises 
from the proper definition of the charge density 

N TV 

a=l 

where the constant TV/2 has been subtracted from the charge density operator in order to make it 
odd under the charge conjugation transformation. As a consequence, when TV is even, p(x) admits 
zero eigenvalues and the ground state does not support any electric flux, while when TV is odd the 
ground state exhibits a staggered configuration of the charge density and electromagnetic fluxes. 

In the continuum the Coleman theorem [ p2[ prevents the formation of either an isoscalar chiral 
condensate < ipip > or an isovector chiral condensate < i[>T a ip > — where T° is an SU (TV) generator 
— for every model with an internal SU (TV)-flavor symmetry. This feature is reproduced on the lattice 
also for this class of models Jl6| . 

4.1 The continuum AT- flavor Schwinger models 

The continuum SU (TV) -flavor Schwinger models are defined by the action 

r M _ 1 

S = / d 2 xC£ V^7^ + 7/^>« - Til V") ( 4 ' 236 ) 
J o=l c 

where the TV fermions have been introduced in a completely symmetric way. Although the theory 



described by ( 4.236 ) strictly parallels what has been shown in the previous section for the SU(2) 
model, we shall now report a detailed analysis both for the sake of clarity and to show that some 
difference appears between TV even and odd. 

The Dirac fields are an TV-plet, i.e. transform according to the fundamental representation of 
the flavor group while the electromagnetic field is an SU(AT) singlet. The flavor symmetry of the 
theory cannot be spontaneously broken for the same reasons as in the SU(2) case. The particles of 
the theory belong to SU (TV) multiplets. The action is invariant under the symmetry 

SU L (Af) <g> SU R (AT) <8> U v (l) ® U A (1) (4.237) 



31 



The symmetry generators act as follows 

(4.238) 
(4.239) 
(4.240) 
(4.241) 

where T a are the generators of the SU(AT) group, 9 a , 9(x) and a are real coefficients and 

Pl = 1(1-k) , P* = ~(l + 75) ■ (4.242) 



SU L {Af) 


1pa(x) ~ 




PL )afc ^6 (») , ip a {x) 


-^- b (x){e-^ Tap «) ba 


SU R {N) 


4>a{x) — 


-> (e^ 


Pa )ab 1pb{x) , "0a (x) 


— (e^ TapL ) ba 


U V (1) 






1 )ah ^b(x) , -0l(x) ~ 




Ua(1) 


1pa{x) ~ 




)o6 06(x) , 





At the classical level the above symmetries lead to conservation laws for the isovector, vector and 
axial currents 



f a {x) R 
jH(x)l 
f{x) 
J^x) 



^ a (x)^P R (T a ) ab Mx) 

^a(x)^P L (T a ) ab Mx) 
tPa(xh' 1 'i-a b 1p b (x) , 
tpa(xh t "l 5 'i-a b tp b (x) . 



(4.243) 
(4.244) 
(4.245) 
(4.246) 



At the quantum level the vector and axial currents cannot be simultaneously conserved. If the 
regularization is gauge invariant, so that the vector current is conserved, then the axial current 
acquires the anomaly which breaks the symmetry 1/4(1) |E3 



( x ) = N -^e^F^(x) 
Zn 



(4.247) 



The isoscalar < ipx/j > and isovector < ipT a ifi > chiral condensates are zero due to the Coleman 
theorem |^2| , in fact they would break not only the Ua{1) symmetry of the action but also the 
continuum internal symmetry SUi,(Af) <8> SU R (Af) down to SUv(Af). There is an order parameter 
just for the breaking of the J7a(1) symmetry J42, Wffl, the operator 



Under a discrete chiral rotation 

0L " 



75 ipi 



-tp L 



CAT) ^_ f e \Mf, / A/ '„ ^AA 



4>R -> JblpR = Ipi 



the operator (4.248) of course transforms as 



^r'>^ (-1)^<^ 



-AM) T (l) (i 



0l m 



> 



(4.248) 



(4.249) 



(4.250) 



The operator (4.248) is even under (4.249) when Af is even and this implies that notwithstanding 
the fact that the continuo us chi ral rotations 1/4(1) are broken by the non-zero VEV (4.248), the 
discret e chiral symmetry ( 1.249 ) is unb roken. When Af is odd also the discrete chiral symmetry 
(|4.249|) is broken by the non-zero VEV (|4.248|) . 

The usual abelian bosonization procedure may again be applied provided that Af Bose fields are 
introduced SI [si El 



l^Wa := -LeTdvQa ,o=l,. 

\/7T 



.Af 



The electric charge density and the action read 

AT 



Jo 



(4.251) 



(4.252) 



a=l 
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N 



N 



S= [ d 2 x{\ £ d M * a d>* a + $a) 2 ) (4.253) 

•* a=l a=l 

The mass matrix is determined by the last term in Eg. ( 4.253] ) and must be diagonalized. The field 
degrees of freedom span the vector space on which the mass matrix is defined. The action must be 
expressed in terms of an orthonormal basis of field vectors, in order to have a properly normalised 
kinetic energy term. The original $ a in Eq.( 4.253) are orthonormal basis vectors, but they are not 

2 

eigenvectors of the mass matrix. The mass matrix has one non-zero eigenvalue jf^ with associated 
eigenvector ^= Y^a=i ® a an d an the other eigenvalues are zero. The remaining eigenvectors can be 
made orthonormal by the following change of variables 



where the orthogonal matrices O b arc p4 



Ol 



Ot 



o: 



1 



(M, 



i 



VWW-l) 

i 



,i) 
(1,1, 



vW-W-2) 



■•,-AA+l) , 
(l,l,...,-AA+2,0) 



(4.254) 

(4.255) 
(4.256) 
(4.257) 



^=(1,-1,0,. ..,0) 



In terms of these new fields $ a the action ( 4.253| ) reads 



S= / d 2 x{-Y^d^ a d^ a 

J 1 a=l 



(4.258) 



(4.259) 



where /x 2 = The action ( 4.259 ) describes M non interacting fields, one massive and M — I 

massless. The multiflavor Schwinger model can also be studied in the framework of non abelian 
bosonization f48[| , where the relationship between isovector currents and bosonic excitations appears 
in a more symmetrical form [^9| . 

The Hamiltonian, gauge constraint and non-vanishing (anti-)commutators of the continuum 
TV-flavor Schwinger models are 

r e 2 M 

H= dx[-E 2 (x)+J2 i/>t(x)a(W x + eA(x))il> a (x)] (4.260) 

J 1 a=l 

d x E(x) + Eli ^aWMx) ~ (4.261) 
[A(x), E(y)} = iS(x - y) , {^ a (x), 4(y)} = S ab S(x - y) (4.262) 



4.2 The lattice AT- flavor Schwinger models 



On the lattice the Hamiltonian, constraint and (anti-) commutators reducing to ( 4. 26C . 4. 261 , 4. 262 ) 
in the naive continuum limit are 



2 N 

eia 



E x — E x _i 



[A X ,E, 



Hs = ^J2E 2 x - |- £ $>l, x+1 e^a,* - i>l,e- iA ^ a<x+1 ) 



N M 



2a 



A'' 



a=l 



(4.263) 
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The fermion fields are defined on the sites, x = 1,...,N, gauge and the electric fields, A x and E x , 
on the links N is an even integer and, when N is finite it is convenient to impose periodic 

boundary conditions. When N is finite, the continuum limit is the A/"-flavor Schwinger model on a 
circle [Q. The coefficient t of the hopping term in (1.262) plays the role of the lattice light speed. 
In the naive continuum limit, = e c and t = 1. 

The lattice A/"-flavor Schwinger model is equivalent to a one dimensional quantum Coulomb gas 
on the lattice with AT kinds of particles. To see this one can fix the gauge, A x — A (Coulomb 
gauge). Eliminating the non-constant electric field and using the gauge constraint, one obtains the 
effective Hamiltonian 



H,, + H.„ 



2N 



E 2 + 2j^J2p(x)V(x-y)p(y) 



■r.y 



it 
2a 



AT 



J2^l*+i eiA ^* ~ ^U e " M ^a,*+l) 



(4.264) 



where p(x) is given in Eq. (4.235) and the Coulomb potential V(x — y) is given in Eq.(3.184). 
The constant electric field is normalized so that [A, E] — i . The constant modes of the gauge 
field decouple in the thermodynamic limit N — > oo. In the thermodynamic limit the Schwinger 
Hamiltonian (4.264), rescaled by the factor e^a/2, reads 

H = H Q + eH h 

with 



H 



E 



{x-yf 



N 

x>y 

—i(R — L) 



(x - y)]p(x)p(y) 



(4.265) 

(4.266) 
(4.267) 



and e = t/e 2 L a 2 . In Eq.(4.267) the right R and left L hopping operators are defined {L = EJ) as 



N 



R 



N AT N AT 



x—l x—1 a—1 x—1 a—1 

On a periodic chain the commutation relation 

[R, L] = 

is satisfied. 



(4.268) 



(4.269) 



When Af is even the ground state of the Hamiltonian (4.266) is the state \g.s. > with p(x) = 
on every site, i.e. with every site half-filled 



Af, 



^^lx^ax\g-s. >= —\g.s. > 



(4.270) 



a=l 



It is easy to understand that p(x) = on every site in the ground state by observing that the 
Coulomb Hamiltonian (4.266) is a non-negative operator and that the states with zero charge 
density are zero eigenvalues of ( 4.266[ ). \g.s. > is an highly degenerate state; in fact at each site x 
the quantum configuration is 



nv-Lio 



> 



(4.271) 



o=l 



— • i.e. it takes on any orientation of 



The state (1.271) is antisymmetric in the indices a = 1, 

the vector in the representation of the flavor symmetry group SU(Af) with Young tableau given in 
fig.fl). The energy of \g.s. > is of order 1, since it is non zero only at the second order in the strong 
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Figure 4: The representation of SU {N) at each site when Af is even 



coupling expansion. 

When AT is odd the ground states of the Hamiltonian ( 4.266 ) are characterized by the staggered 
charge distribution 

p{x) = ±h-\T (4.272) 



since (4.272) minimizes the Coulomb Hamiltonian (4.266); one can have p{x) = +1/2 on the even 
sublattice and p(x) — —1/2 on the odd sublattce or viceversa. The electric fields generated by the 



charge distribution (4.272) are 



Since now 



H \g.s. >= Y^\g-s- > 



(4.273) 



(4.274) 



the ground state energy is of order e\. The states \g.s. > are highly degenerate since they can take 
up any orientation in the vector space which carries the representation of the SU(Af) group with 
the Young tableaux given in fig.(|5|). 

Either when AT is even or when Af is odd the ground state degeneracy is resolved at the second 
order in the strong coupling expansion. First order perturbations to the vacuum energy vanish. 
The vacuum energy at order e 2 reads 



E {2) -< 



n 



H 



(4.275) 



where the expectation values are defined on the degenerate subspace of ground states and II is 
a projection operator projecting orthogonal to the states of the degenerate subspace. Due to the 
commutation relation 



[H 0) H h ] = ~~irj~~ Hh - 2 £ [V(x - y) - V(x - y - 1)] (L y + R v )p{x) 



N 



(4.276) 



Eq.(4.275) can be rewritten as 



E, 



(2) _ 



~2<RL> 



(4.277) 
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Figure 5: The representation of SU(Af) at each site of the even sublattice and odd sublattice when 
Af is odd 



On the ground state the combination RL can be written in terms of the Heisenberg Hamiltonian 
of a generalized SU(Af) antiferromagnet. By introducing as in the previous section the Schwinger 
spin operators 

S x = (4-278) 
where T a are now the generators of the SU(Af) group, the SU(Af) Heisenberg Hamiltonian reads 



N Af 1 1 N 

Hj = S x ■ S x+1 - J - + ^P{x)p{x + !)) = --£ L X R 



x=l 



(4.279) 



When Af is even, on the degenerate ground states one has 

N 



N 



< Hj >=< J2(S X ■ S x+1 - — ) >=< LxRx > 

x—1 x—1 

while when Af is odd one has 



iV 



<Hj >=<J2(Sx-S x+1 



Af 2 + 1, 
87V ' 



1 N 

>=< —-x ^ L X R X > 



(4.280) 



(4.281) 



Taking into account that the products of L x and R y at different points have vanishing expectation 
values on the ground states and using Eq.( 4.280| ) or Eq.( 4.281 ), Eq.( 4.277 ) reads 



EP =4<H.,> 



(4.282) 



The problem of determining the correct ground state, on which to perform the perturbative ex- 
pansion, is then reduced again to the diagonalization of the SU(Af) Heisenberg spin-1/2 Hamil- 
tonian (4.27E). As we already pointed out in section (2.4), generalized SU(Af) antiferromagnetic 
chains have not been yet analysed in the literature in such a detailed way as the SU(2) chains. 
Consequentely, the study of the lattice SU (Af) flavor lattice Schwinger models become extremely 
complicated for Af > 2. Nonetheless, the computational scheme that we developed for the U(l)- 
and 5't/(2)-fiavor models |l4|, [l5|, [ll|, should work for a generic SU(Af)-ft&vor model. 

The ground state of the gauge models is very different depending on if Af is even or odd. When 
Af is even, the ground state \G.S. > of the spin Hamiltonian (4.279) is non-degenerate and transla- 
tionally invariant, and since it is the ground state of the gauge model in the infinite coupling limit, 



3G 



there is no spontaneous breaking of the chiral symmetry for any SU (27V)-flavor lattice Schwinger 



model. In contrast, when Af is odd, the ground state \G.S. > of the spin Hamiltonian (4.279) is 
degenerate of order two and is not translationally invariant and consequently any SU(2Af+ l)-flavor 
lattice Schwinger model exhibits spontaneous symmetry breaking of the discrete axial symmetry. 
By translating of one lattice spacing \G.S. > one gets the other one. The A/"- flavor lattice Schwinger 
models excitations are also generated from \G.S. > by two very different mechanisms, that we al- 
ready described for the two-flavor model in section (3.2). There are excitations involving only flavor 
changes of the fermions without changing the charge density p(x) which corrispond to spin flips 
in the SU(2) invariant model. These excitations are massless. At variance massive excitations in- 
volve fermion transport besides flavor changes and are created by applying to \G.S. > the latticized 
currents of the Schwinger models which vary the on site value of p(x). 

Very different is the case of the massive multiflavor Schwinger models J|5| . When Af is odd, the 
presence of a non-zero fermionic mass m removes the degeneracy and selects one of the two \G.S. > 
as the non degenerate ground state. When Af is even the ground state remains translationally 
invariant in the strong coupling limit e| 3> m. In the weak coupling limit m 3> e 2 L the discrete 
chiral symmetry is broken for every Af. 



5 Concluding remarks 

In these lectures we analysed the correspondence between the multiflavor strongly coupled lattice 
Schwinger models and the antiferromagnetic Heisenberg Hamiltonians to investigate the spectrum 
of the gauge models. 

Using the analysis of the excitations of the finite size chains, we showed the equality of the 
quantum numbers of the states of the Heisenberg model and the low-lying excitations of the two- 
flavor Schwinger model. We provided also the spectrum of the massive excitations of the gauge 
model; in order to extract numerical values for the masses, we explicitly computed the pertinent 
spin-spin correlators of the Heisenberg chain. 

The massless and the massive excitations of the gauge model are created from the spin chain 
ground state with two very different mechanisms: massless excitations involve only spin flipping 
while massive excitations are created by fermion transport besides spin flipping and do not belong to 
the spin chain spectrum. As in the continuum theory, due to the Coleman theorem fl22]j , the massless 
excitations are not Goldstone bosons, but may be regarded as the gapless quantum excitations of 
the spin-1/2 antiferromagnetic Heisenberg chain J56fl , 

In computing the chiral condensate one can show |l6| that, also in the lattice theory, the 
expectation value of the umklapp operator is different from zero, while both < tjjtp > and < ipa°"ip > 
are zero to every order in the strong coupling expansion. This implies that both on the lattice and 
the continuum the SU(2) flavor symmetry is preserved whereas the J7a(1) axial symmetry is broken. 
The umklapp operator jl^, [l6| is the order parameter for this symmetry, but being quadri-linear 
in the fcrmi fields, is invariant, in the continuum, under chiral rotation of n/2 and on the lattice 
under the corresponding discrete axial symmetry (translation by one lattice site) . This shows that 
the discrete axial symmetry is not broken in both cases. Our lattice computation enhance this 
result since the ground state of the strongly coupled two-flavor Schwinger model is translationally 
invariant. 

The pattern of symmetry breaking of the continuum is exactly reproduced even if the Coleman 
theorem does not apply on the lattice and the anomalous symmetry breaking is impossible due to 
the Nielsen-Ninomiya [^7], theorem. At variance with the strongly coupled one-flavor lattice 
Schwinger model, the anomaly is not realized in the lattice theory via the spontaneous breaking 
of a residual chiral symmetry [JL4j, but, rather, by explicit breaking of the chiral symmetry due 



to staggered fermions. The non- vanishing of the VEV of the umklapp operator [ fl5[ 16 may be 
regarded as the only relic, in the strongly coupled lattice theory, of the anomaly of the continuum 
two-flavor Schwinger model. It is due to the coupling induced by the gauge field, between the right 
and left-movers on the lattice. 
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When the fermion mass m is different from zero, some further difference arises between Af odd 
and Af even. When Af is odd, the mass term induces a translational non invariant ground state, 
generating a spontaneous chiral symmetry breaking. When Af is even, the ground state remains 
translationally invariant in the strong coupling limit, i.e. e 2 3> m. In the weak coupling limit, 
m 3> e 2 , the discrete chiral symmetry is spontaneously broken for every Af. For Af = 2, the 
soliton-antisoliton excitations pq acquire a mass. 
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